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In this paper, we present a theoretical description to the quantum coherence and decoherence 
phenomena of energy transfer in photosynthesis observed in a recent experiment [see Science 316, 
1462 (2007)]. As a successive two-color laser pulses with selected frequencies cast on a sample 
of the photosynthetic purple bacterium Rb. sphaeroides two resonant excitations of electrons in 
chromophores can be generated. However, this effective two-level subsystem will interact with 
its protein environment and decoherence is inevitable. We describe this subsystem coupled with 
its environment as a dynamical spin-boson model. The non-Markovian decoherence dynamics is 
described using a quasi-adiabatic propagator path integral (QUAPI) approach. With the photon- 
induced effective time-dependent level splitting energy and level flip coupling coefficient between 
the two excited states and the environment-induced non-Markovian decoherence dynamics, our 
theoretical result is in good agreement with the experimental data. 
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I. INTRODUCTION 

In the past decades, increasing attention has been paid 
to understand the intrinsic mechanisms of the efficient 
energy transfer in Hght harvesting complexes. In the 
previous investigations, the energy transfer was often 
described by a semiclassical model that invokes 'hop- 
ping' excited-state populations along discrete energy lev- 
els. Explicitly, the electrons in pigment molecules, chro- 
mophores, are excited by the incident light. Then the 
electronic excitation moves downhill from energy level 
to energy level through the chromophores before being 
trapped in the reaction center. However, it has not been 
conceived that the high efficient energy transfers between 
chromophores is realized through the electronic excited 
state hopping. 

To understand the energy transfer among excited 
states of each chromophore in the Rhodobacter (Rb.) 
sphaeroides reaction centers, Fleming's group proposed 
recently that the plants and the photosynthetic bacteria 
may utilize a 'clever' quantum algorithm for the energy 
transfer [l|, Q . Unlike the semiclassical model, Fleming 
and co-workers suggested that a successive excitations 
make a coherence superposition between electronic ex- 
cited states of the chromophores. The energy transfer 
is considered to be carried by the coherence of the su- 
perposition state through an easiest route analogous to 
Grover's algorithm for the fastest possible search of an 
unsorted information database 0]. However, Mohsehi et 
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al. [3] argued that the standard Grover's algorithm may 
not be applicable to the energy transfer mechanism in 
photosynthesis. Nevertheless, even if the plants and bac- 
teria work unlike the Grover's algorithm, the coherent dy- 
namics evolution of the electronic excited states in chro- 
mophores could still play an important role in the chro- 
mophoric energy transfer [1, @ . In this paper we shall de- 
rive a dynamical spin-boson model to describe how long 
the coherent superposition state can be persisted in the 
photosynthetic purple bacterium Rb. sphaeroides. 



II. DYNAMICAL SPIN-BOSON MODEL 

Recent experiments by Fleming group showed that the 
electronic coherence between two excited states may last 
for 400 fs or longer [l|, [2]. Instead of investigating the 
coherence among all excited states of each chromophores 
in the Rb. sphaeroides reaction centers 0,[1|, they inves- 
tigated the coherent superposition state of two electronic 
excited states \H) and \B) in chromophores created ex- 
perimentally from two successive laser pulses. The laser 
pulses (with wavelength 750 nm and 800 nm) derives the 
system which is initially in the ground state \g) into a co- 
herent superposition state \'^) ~ ct\H) + l3\B). However, 
the system of the two excited states inevitably interacts 
with its environment, which results in the decoherence 
of the coherent superposition state j^"). To obtain the 
persistent time of coherence, i.e., the decoherence time 
of I^P), Fleming et al. use the third incident pulse (with 
wavelength 750 nm) to produce a photon echo. From 
the echo signals they extract the decoherence time of the 
coherent state jvl/) 

To be specific, in the two-color photon echo experi- 
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ment on bacterial reaction centers (RC) the RC from 
the photosynthetic purple bacterium Rb. sphaeroides in- 
cludes a bacteriochlorophyll dimer, called the special pair 
(p), in the center, an accessory bacteriochlorophyll flank- 
ing p on each side (BChl), and a bactcriopheophytin 
(BPhy) next to each BChl. The RC of Rb. sphaeroides 
has several absorption peaks which are made by the chro- 
mophores. In Ref. the absorption spectrum of the 
p-oxidized RC at 77 K shows the H band at 750 nm and 
the B band at 800 nm (where H and B are used to de- 
note excitonic states which are dominatingly produced 
from monomeric BPhy and accessory BChl in the RC, 
respectively). At the first step, they used a successive 
two laser pulses with different colors and tuned for res- 
onant excitation of H transition at 750 nm (at time r{) 
and the B transition at 800 nm (at time T2) cast on the 
sample. The two pulses produce a coherent superposi- 
tion state 1 5*) between the electronic excited states \H) 
and \B). Then, in order to measure the decoherence time 
of the coherent superposition of these two excited states, 
the third laser pulse is cast on the sample after a time 
t2 from the second pulse, which generates a photon echo. 
When the times ti = — t[ and t2 = — are dif- 
ferent one can detect the different integrated intensity of 
the echo signals in the phase matched direction. If the 
time ti is fixed (the fixed time ti = 30 fs in Ref. [l| ), the 
integrated echo signals as a function of t2 represents the 
decoherence of the coherent superposition between \H) 
and \B). 

The evolutions of the integrated echo signals are plot- 
ted in the Fig. 3 of Ref. [l| . From the experimental data 
one can see that: (1) The coherence is resonant with dif- 
ferent frequencies in the first 400 fs. The first, the second 
and the third periods are about 100 fs, 120 fs, and 130 
fs respectively; (2) The coherence between the excited 
states \H) and \B) persists for more than 400 fs; (3) The 
third peak in the evolution of the coherence decay more 
fast than the first and second peaks. In the following, 
we shall attempt to understand theoretically the dynam- 
ical process of excited electron states in chromophores 
and to explain the above experimental result. It is in- 
teresting to see, as we shall show later, that the prob- 
lem can be described by a dynamical spin-boson model 
where the two excited states are time-dependently cou- 
pled one another through the laser pulses, as an effect of 
the photon-induced dynamics. To be explicit, we model 
the BChl and BPhy molecules by the Hamiltonian in the 
Condon approximation as [H, [Tj, [lB| 

He =eo \g) {g\ + (H \H) {H\ + (B \B) {B\ 

+ eHB\HB){HB\ + Jo {\H) {B\ + \B) 

+ H-E{t){\g){j\ + \j){g\). (1) 

3=H,B 

Here, eo, tn, and chb are the energies of the ground 
state I (7), the two excited states \H) and \B), and the 
doubly excited state \HB). Jq is the electronic couphng 
between the two excited states, flj (j = H, B) is the cor- 



responding electronic dipole, and E (t) is the external 
electronic filed of the two successive laser pulses, namely, 
E{t) = E^^e-'^iit+tif+iujHit+ti) j^^^^~r2t^+iujBt ^j^^j-g 

ri,2 are the decay constants of the laser pulses as they 
pass though the sample. 

For simplicity, we assume that the dipole moments for 
the H and B states are the same: fin = fln = We can 
decouple the states \g), \HB) from \B) by making 
the following canonical transformation to Eq. ([1]): 



S = Kit)[ai\H) (5I-I.9) {H\) 
+ P{\B) {g\-\9) {B\)], 



(2) 
(3) 



where K{t) — jl- E{t), a ~ {tB ~ eo ^ Jo)/^, P = {^h — 
eo — Jo)/^, and ft = {en ~ eo)(es — eo) — Jo- Under the 
condition K(t)a ~ K{t)P <ti 1, we obtain 

We = {e^, 1.9) (5I + eHB\HB){HB\} + {e^ \H) {H\ 

J,sm{B\ + \B){H\)}, (4) 



\B) {B\ 



where, eg 



= eo - K'^{t){a + /3), c'b ^ eH + K^{t)a, 
e^ = es -h K'^{t)l3 and J^s = Jo + K^{t){a + /3)/2. As 
one can see, the ground state \g) and the doubly excited 
state \HB) are now decoupled from the single excited 
states \H) and \B). Thus the decoherence dynamics of 
the coherent superposition state is fully determined 
by the effective two-level Hamiltonian, the second curly 
bracket in Eq. ([4]) which can be rewritten as 



A(<) 



(5) 



Here, Ui {i = x, z) are the Pauli matrix, e (<) = e^f — 
(B + K^(<) [a - /3), A (t) = 2Jo + K^{t){a + /?). It shows 
that the energy splitting of the two excited states and the 
coupling between them are shifted by the pulse-induced 
time-dependent dipole-dipole interaction [~ (i)] , as an 
effect of photon- induced dynamics. 

Furthermore, this two-level system, as a part of chro- 
mophores, is inevitable to interact with its protein envi- 
ronment through the thermal vibrations. This thermal 
reservoir can always be modeled with a set of harmonic 
oscillators. Then the total system of the two-level system 
coupled to its environment can be written as a dynamical 
spin-boson model. 



H = iJo(t) + ^ c^{h\ + h) + y^ hLO,b%, (6) 



in which the last two terms are the interaction between 
the two-level system and the bath and the Hamiltonian 
of the thermal bath itself, respectively. The parameter c; 
is the coupling of the system to the bath and b\ [hi) are 
the creation (annihilation) operators of the i-th thermal 
mode. The decoherence dynamics of the above two-level 
system are mainly induced by the back-actions of the 
thermal bath as well as the decay of the time-dependent 
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coupling induced by laser pulses. The solution of the 
problem is now completely determined by the parameters 
e{t),A{t) and the spectral density of the thermal bath[llj| 
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6{lU - LU,) 



u/u 



(7) 

Here ujc is the high-frequency cut-off of the bath modes. 
Different s correspond to the super-Ohmic (s > 1), the 
Ohmic (s = 1), and the sub-Ohmic (0 < s < 1) baths, 
and is a dimensionless dissipative parameter describing 
the coupling between the two-level system and the bath. 



III. NUMERICAL METHOD 
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The spin-boson model has been investigated by many 
methods though it has not be exactly solved and some 
approximation, such as the Markov approximation, must 
be used Here, we shall apply the quasiadiabatic 

propagator path integral (QUAPI) technique [l^ to ex- 
plore the decoherence dynamics of this system, where 
non-Markovian processes are involved. Indeed, as it has 
been pointed out f]\ the evolution of this time-dependent 
system should be highly non-Markovian. To see how the 
non-Markovian dynamics may play an important role, 
it is useful to estimate the correlation time of the ther- 
mal bath, which can be obtained from the bath response 
function 



1 r 

C{t) = - / 



dijjJ {ijj) 



coth 



cos Lut — i sin ujt 



(8) 

where /3 = l/fcsT with Boltzmann's constant ks and 
the temperature T. When the real and imaginary parts 
of C{t) behave as a delta function S (t) and its derivative 
5' (t) , respectively, the dynamics of the reduced density 
matrix is Markovian. Otherwise, non-Markovian dynam- 
ics occurs. The broader the Re[C (t)] and Im[C {t)] are, 
the longer the correlation time will be, and the more se- 
rious the practical dynamics is distorted by the Markov 
approximation. Similar to Ref. [l7|, we calculate the cor- 
relation times of the bath. The result is shown in Fig. 1. 
As we can see the correlation time of the bath is about 
Tc ~ 15 fs which actually depends on the frequency cutoff 
LUc in the bath. 

The decoherence dynamics of the coherent superposi- 
tion states 1 4*) is characterized by the time evolution of 
the reduced density matrix obtained after tracing out the 
bath degrees of freedom, i.e.. 



p(i) = Trbath[e-'^*/''ptot (0) e'^*/''] 



(9) 



Following the experiment, the interaction between sys- 
tem and bath is turned on at t = T2 when the sec- 
ond laser pulse is applied. Thus, the density matrix 
of the total system before the time t = is a direct 
product of the system and bath components, namely, 
Ptot (T2) = p (tj) (g> /Obath (T2) , where p (r^) and pbath (T2) 



FIG. 1: (Color online) Real and imaginary parts of the re- 
sponse functions {C{t)) for Ohmic bath. Here, we take 
= 0.75, hwc = 2000 cm-i and T" = 77 K. 



are the "initial" states of the system and the bath at 
t = Tj. If we set -Boi = £-02 = Eq, the initial resonant 
two excited state can be written as \'^) = (1^) + 1-^))- 
After the canonical transformation, the initial state re- 



mains almost the same: |*') = e'^ 1*) w {\H) + \B)) 

under the condition Ky^ o? -I- fP- <C 1, which is satisfied 
for the parameters we taken in the following numerical 
calculations. We also set the bath initially at the thermal 
equilibrium, namely, pbath ('''2) — e^^^^ j^xieT^^^^^. 

The reduced density matrix p(t) can be evaluated by 
using the well established iterative tensor multiplication 
(ITM) algorithm derived from the QUAPI. This algo- 
rithm is numerically exact and successfully tested and 
adopted in various problems of open quantum systems 
(TgI [Tsj . For details of the scheme, we refer to previ- 
ous works The QUAPI asks for the system Hamil- 
tonian splitting into two parts iJg and H^nv, where 
Henv = He-b + Hb. In Order to make the calculations 
converge we use the time step At = 5 fs which is shorter 
than the correlation time of the bath and the charac- 
teristic time of the two-level subsystem. In order to in- 
clude all non-Markovian effect of the bath in the ITM 
scheme, one should choose A/ci„ax so that A/cmaxAi is 
not much shorter than the correlation time of the bath. 
Here, A/c„iax is roughly equal to the number of time steps 
needed to span the half-width of the response function 
C{t — t') [I4I. Then taking Afcmax = 3 should be large 
enough in our calculations. 
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FIG. 2: (Color online) The evolutions of the off-diagonal co- 
herent terms of the reduced density matrix for the two-level 
subsystem with the Ohmic bath at temperature (a) 77 K and 
(b) 180 K. 



IV. RESULTS AND DISCUSSIONS 

The decoherence of the two exciton states is reflected 
through the decays of the off-diagonal reduced density 
matrix element. In Fig. 2a and b we plot the evolution of 
the off-diagonal reduced density matrix element at tem- 
perature r = 77 K and T = 180 K, respectively, where 
the environment is assumed as an Ohmic bath (s = 1), 
and the parameters en — 12108 cm^^, = 12000 cm~^, 
eo = 10570 cm"\ Jo = 20 cm-\ kq = fl-Eo = 210 cm~^ 
and Ti = = r = 3 X lO^"* s'^. We take the Kondo 
parameter — 0.75 and the frequency cutoff hujc = 2000 
cm"-'^ [ill for the thermal bath. According to Ref. we 
take ti = 30 fs, and ^2 — 50 fs. The black square boxes in 
Fig. 2 are the experimental data from [l| , the linked blue 
and red dots are our numerical results. The plots show 
that our theoretical description is in good agreement with 
the experimental results. As we see the coherence decay 
is much faster after the second peak and the oscillation 
periods of the evolutions increase with the time. The os- 
cillation dephasing behaviors indeed come from the time- 
dependent level splitting e (t) and coupling A (t) induced 
by laser pulses as well as the non-Markovian processes 
due to the interaction between the two-level system and 
the thermal bath. 

Further calculations show that using the sub-Ohmic 
bath with £^1/2 = 0.07 and super-Ohmic bath with ^2 — 
50, we can obtain the similar results, as show in Fig.jHK- 
However, if we use the same value of the dimensionless 
dissipative parameter for the Ohmic, sub-Ohmic and 



FIG. 3: (Color online) The evolutions of the off-diagonal co- 
herent terms of the reduced density matrix for the two-level 
subsystem in sub-Ohmic and super-Ohmic Baths. Here, we 
set 5i/2 = 0.07 and ^2 = 50 at T = 77 K. Other parameters 
are the same as that in Fig. 2 



super-Ohmic spectral densities, the corresponding deco- 
herence behaviors are completely different. This indi- 
cates that the decoherence dynamics is non-Markovian. 
To show the non-Markovian effect, we also calculate the 
Markov approximation with the same parameters in the 
Ohmic case. The detailed numerical results are plotted in 
Fig.[3|D, where the difference between the non-Markovian 
dynamics and the corresponding Markov approximation 
is obvious. 

We should also point out that the bare level splitting 
e/f — es is not the experimentally observed values. The 
experimental values 600 cm~^) of the level splitting 
contains the pulse-induced time-dependent effects plus 
the back-action effect from the environment in the non- 
Markovian regime so that the input value of en — £b 
is different from the splitting of experimental linear ab- 
sorbance peaks. A practical calculation of the true site 
energy difference that contains all the the pulse-induced 
time-dependent effects plus the back-action effect from 
the environment may be possible [l^] but it is beyond 
the scope of the present investigation. Here the input 
value of the bare level splitting is based on how better 
to fit the experimental data. Our numerical calculations 
show that the profiles of the evolution curves for the off- 
diagonal reduced density matrix elements are sensitive to 
the changes of the bare level splitting £// — es, the bare 
coupling Jo and the parameter kq but insensitive to the 
changes of ch and es with €h — £b ~ HO cm~^. In Fig. 21 
we plot the time evolution of off-diagonal reduced density 
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FIG. 4: (Color online) The evolutions of the off-diagonal co- 
herent terms of the reduced density matrix in Ohmic bath, 
with different bare level splitting en — es = 108, 118, 158 cm~^ 
at T — 77 K. Other parameters are the same as that in Fig. 
2 



FIG. 5: (Color online) The evolutions of the off-diagonal co- 
herent terms of the reduced density matrix for the two-level 
subsystem in Ohmic bath at T = 77 K, with different bare 
electronic coupling Jo = 20,5, and 0.01 cm~^. Other param- 
eters are the same as that in Fig. 2. 



matrix element with different en — ^b- As we can see, 
the change of the bare level splitting en ~ produces a 
quite different pattern of the decoherence dynamics. A 
smaller value of en — £b will result in a shorter oscillation 
period of the evolution in time. 

Also, the bare electronic coupling Jq between the exci- 
ton states \B) and \H) is not the electronic coupling J be- 
tween the BChl and BPhy molecules that used in [l| . Our 
numerical results suggest that Jq is much smaller (< 30 
cm^^) than the pulse-induced dipole-dipole interaction 
but is not negligible. Change of Jq leads to very different 
decoherence behavior as shown in Fig. [5l which is actu- 
ally also convinced in the supplemental material of [^. 
The resulting effective electronic coupling Jeff oscillates 
in time with the maximum amplitude being less than 
400 cm~^. As one can see, our numerical result shows 
that the oscillation of the off-diagonal reduced density 
matrix element is due to the effective time-dependence 
of the level splitting and flip coupling between the two 
excitonic states induced by the laser pulses. This differs 
from the interpretation given in 1] where it has been 
pointed out that reproducing the long dephasing time 
observed would require a value of J > 320 cm~^, but a 
large J implies that the BChl and BPhy excitations are 
almost degenerate, which may not be consistent with the 
observed spectra. Therefore they argued that the vibra- 
tional modes begin to contribute in the energy transfer. 
Their simulation shows that a vibrational mode (with a 
frequency to = 250 cm~^) coupling to the excitonic states 



can fit the long dephasing time very well. Here we show 
that the large effective coupling between the two excited 
states induced by the laser pulses combining with a very 
smaller but not negligible bare coupling Jq can reproduce 
the long dephasing time observed, without including the 
contribution arisen from the vibrational mode. 

On the other hand, the level splitting and flip cou- 
pling we derived explicitly depend on the parameter 
kq — fl ■ Eq- The dipole moment fj, = ^/D (debye) 
where the dipole str eng th D can be estimated by: D ~ 
0.0196nemax<5/A„iax UM- Here n = 1.359 is the refractive 
index, emax and S are the half-width and the peak value at 
Amax = 750 nm and 800 nm for H and B excited states, 
respectively. Fig. 1 of fl'] shows that e^^^^ — 0.5eJ^Q^ 
and 5p ~ 1.1 (5b- Taking Z3b = 40 debye^ (^corre- 
sponding to fl — 6.3 debye for 800 nm BChl [22]), we 
have Z^H — 23.5 debye^. Averaging the dipole strength 
D = (40-1-23. 5)/2 for BChl and BPhy, we obtain fj. ~ 5.63 
debye = 1.88 x 10"^^ C cm. The power intensity used in 
Fleming's experiment [J is P = 1.3 x 10~^ J/cm^ while 
the pulse duration At = 40 fs. Therefore the pulse inten- 
sity is Ip — P/ At — 3.25 X 10^ W/cm^ which corresponds 
to a pulse field amplitude Eq = V^Erm — 2.21 x 10^ 
V/cm. Thus, the realistic dipole-field couphng strength 
kq = /i • i?o — 209 cm^^ . This is very close to the value of 
our theoretical best fitting kq = 210 cm~^ in the numeri- 
cal calculation. We also find that the laser pulse strength 
Eq controls the decay curves. In Fig. [SI we plot the evo- 
lution of the off-diagonal reduced density matrix element 
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periods of the evolution in time, and vice versa. This 
property may be used for further experimental test of 
whether the pulse-induced time-dependent dipole-dipole 
interaction or the additional vibrational mode gives rise 
to the oscillation decay of the coherent cxcitonic states. 

In conclusion, we present in this paper a dynamical 
spin-boson model with the effective time-dependence of 
the level splitting and flip coupling between the two ex- 
cited states, induced by laser pulses, to describe the long 
dephasing time observed recently in a photosynthesis ex- 
periment [ij] . We use the quasi-adiabatic propagator path 
integral approach to analyze in detail the non-Markovian 
decoherence behaviors of the two excitonic states, \H) 
and \B). Due to the photon- induced time-dependence 
of the energy level splitting and flip coupling coefficient 
between the two excited states and the environment- 
induced non-Markovian dynamics, our theoretical result 
is in good agreement with the experimental data. 



FIG. 6: (Color online) The evolutions of the off-diagonal co- 
herent terms of the reduced density matrix for the two-level 
subsystem with different pulse field strength, kq = /2 ■ -Eo = 
170,210 and 250 cm"\ in Ohmic bath at T = 77 K . Other 
parameters are the same as that in Fig. 2 

to demonstrate the KQ-dependence of the dephasing time. 
It shows that increasing Eq will decrease the oscillation 
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